ON SIMILARITY SOLUTIONS TO THE 
MULTIDIMENSIONAL AGGREGATION EQUATION 



HONGJIE DONG 

Abstract. We study similarity solutions to the multidimensional 
aggregation equation u t + drv(uv) = 0, v = — VA * u with general 
power-law kernels A such that VAT(x) = x\x\ a ~ 2 ,a G (2 — d,2). 
We analyze the equation in different regimes of the parameter a. 
In the case when a g [4 — d, 2), we give a characterization of all the 
"first kind" radially symmetric similarity solutions. We prove that 
any such solution is a linear combination of a delta ring and a delta 
mass at the origin. On the other hand, when a e (2 — d, 4 — d), 
we show that there exist multi delta-ring similarity solutions in 
R d . In particular, our results imply that multi delta-ring similarity 
solutions exist in 3D if a is just a little bit below 1. 



1. Introduction 
In this paper, we study the multidimensional aggregation equation 

u t + div(uv ) = 0, v = — VK * u (1.1) 

for x G M. d and t > with the initial data 

u(0,x) = u (x), x e R d . 

Here d is the space dimension, u > is the density function, v is the 
velocity field, K is the interaction potential kernel which typically is 
taken to be radially symmetric, and * denotes the spatial convolution. 
The aggregation equation arises in various models for biological ag- 
gregation and problems in granular media; see, for instance, Mogilner 
and Edelstein-Keshet [16] and Carrillo, McCann, and Villani [8]. The 
equation appears as an alternative to the Debye-Huckel or Keller- 
Segel model. It can be viewed as a continuum model for N particles 
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Xi, . . . , Xm-, which interact via a pairwise interaction potential: 
rlX N 

^ = -J2m k VK(X l -X K ). 

k=l 

Examples of such kernels which are commonly used are K(x) = 1 — e~' x ' 
and K(x) = \x\ a . For nonsmooth kernels K, a typical property of the 
solutions is that they concentrate mass in finite time. On the other 
hand, for smooth kernels, solutions blow up at the infinity (cf. [I]). 

In the past few years, the aggregation equation has attracted many 
attentions due to its diverse applications in biology and physics. The 
problems of the well-posedness in various spaces, the existence of finite- 
time blowups, asymptotic behaviors of solutions of this equation have 
been studied extensively in a series of papers; see, for example, Li 
and Toscani [T5], Bodnar and Velazquez [B], Bertozzi and Laurent [2], 
Laurent [12], Bertozzi, Laurent, and Rosado [I], Carrillo et al. [TJ, the 
author [9], Bertozzi, Garnett, and Laurent j5] and references therein. 
The equation with an additional dissipation term is closely relevant 
to the Keller-Segel model and was investigated by Li and Rodrigo in 
[T3t IT3] . See also [3] for a nice review about this equation. 

We are interested in profiles of first-kind self-similar solutions to 
the equation (11. ip with general power-law interaction kernel K(x) = 
\x\ a /a, a G (2 — d, 2) (when a = 0, naturally we take K(x) = log \x\). 
These solutions are mass conserving and of the form 

u (^ x ) = T^ u o ( 7^77 J • (I- 2 ) 



R(t) d u \R(t) 

For many evolutionary equations, this type of solutions are of special 
interest as they usually give large time behaviors of global solutions 
or asymptotic behaviors of finite-time blowups evolving from general 
initial data. One class of such solutions is the trivial single delta-ring 
solutions, where u = m5\ x \ =p for m,p > 0, i.e., the support of the so- 
lution concentrated on a sphere with shrinking radius. Recently in [T] , 
Bertozzi, Carrillo and Laurent studied nontrivial similarity solutions 
to (II. ip with a = 1. In ID, they constructed a weak solution with 
support on an open interval: 



u{t,x) = — -U 



x 



t \T* - t 

where U is the uniform distribution on (—1,1). While in 2D, they 
found a double delta-ring solution given by 



u = mi6\ x \ =Pl + m 2 S\ x \-_ 



-in 
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for appropriate rrii, pi > 0. In the same paper, using a relation between 
the kernel K and the Laplace operator, they also proved that in any 
odd dimension d > 3 such solutions cannot exist with support on open 
sets. Later in [9] , by observing certain concavity property of the kernel 
in the polar coordinates, we completely characterized all the radially 
symmetric first-kind self-similar solutions in any dimension d > 3 when 
a = 1: all such solutions must have the form (11.21) . where R(t) is a 
linear positive decreasing function and 

u = m S\ x \ =0 + mi<5| x | =pi 

for any mQ,mi > 0, p\ > 0. In particular, this result implies that with 
a = 1 the equation (11.11) does not have any double or multi delta-ring 
solutions when d > 3. We also mention that in a recent interesting 
paper [11], Huang and Bertozzi provided very precise numerical evi- 
dence of the local existence of smooth radially symmetric similarity 
solutions to (II. ip when a = 1 in high space dimensions, which are of 
the second kind and do not conserve mass. Similar numerical simu- 
lations for a G (0, 2) have been carried out by the same authors in a 
recent preprint |10j . in which they also studied the long time behavior 
of global smooth similarity solutions when a > 2. 

The objective of this paper is to investigate similarity solutions to 
(11.11) with general power-law kernel by using analytic methods. To the 
best of our knowledge, except for some numerical results mentioned 
above, this problem has not been studied in detail. We analyze the 
equation in different regimes of the parameter a. Since VK = x\x\ a ~ 2 , 
the velocity field is always contracting. In the case when a G [4 — d, 2), 
we give a characterization of all the first-kind radially symmetric simi- 
larity solutions. We prove that any such solution is a linear combination 
of a delta ring and a delta mass at the origin; see Theorems 12.21 and 
13.11 This result generalizes the aforementioned results in [1] and [9] 
proved for the case a = 1. On the other hand, when a G (2 — d, 4 — d), 
we show that there exist multi delta-ring similarity solutions in IR d ; see 
Theorem 14.11 This, in particular, implies the existence of multi delta- 
ring similarity solutions in 3D if a is just a little bit below 1. Finally, 
we present an example of first-kind similarity solutions with support 
on open sets in the limiting case a = 2 — d, which generalizes the ID 
example in [T] with a = 1 ; see Example 15.11 

All the proofs in this paper are based on delicate analysis of the 
kernel in the polar coordinates; see (I2.3P for the definition of 0. For 
the proofs of Theorems 12.21 and I3.1[ first we note that the arguments 
in [1] and [9] break down in several places for general a. In particular, 
there is no clear relation between K and the Laplace operator, and we 
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do not have the concavity of the kernel in the polar coordinates; cf. 
Remark 12.41 Here we refine an idea in [S] by observing that instead 
of the concavity of 0, the proof of the corresponding result in [9] only 
uses the monotonicity property of 0/r. This enables us to extend the 
argument to the general a G [4 — d, 2). However, we are not able to 
find a unified proof of this monotonicity property. We treat the case 
« = 4 — d or a > 5 — d in Section [2] and the case a G (4 — d, 5 — d) in 
Section [3] separately by using completely different methods. The proof 
of the existence of multi delta-ring solutions (Theorem 14. lj) relies on 
the fact that 0'(O) < 0(1) when a G (2 — d, 4 — d). Once this is verified, 
we reduce the problem to the existence of positive solutions to a system 
of linear algebraic equations. 

Regarding the double and multi delta-ring similarity solutions con- 
structed in this paper and [T], an interesting question would be whether 
these solutions describe the asymptotic behavior of certain finite-time 
blowup solutions. Another relevant problem is the stability of these 
solutions. We shall study these problems in future work. 

2. Characterization of similarity solutions 

Let u be a radially symmetric first-kind similarity solution given by 
(11. 2p . Recall that we take VK = x\x\ a ~ 2 . By the homogeneity, it is 
easily seen that 



Upon using a similar argument as in |T] by looking at the flow map 
driven by v , we know that on the support of Uq, 



with some constant A > 0. Indeed, the function u given by (II. 2p is 
a first-kind similarity solution of (11. ip if and only if for any x on the 
support of Uq, the flow map starting from x is equal to R(t)x at time t, 
which clearly implies (12. ip . Moreover, since R(t) satisfies the ordinary 
differential equation 



where To = (2— a)~ 1 A~ 1 is the first time when all the mass concentrates 
at the origin. 

We denote &(M. d ) to be the set of all probability measures on M d . 




Vq = —X\X 




Vq = — Xx 



(2.1) 



R'(t) = v(t,R(t)) = -\R(t) 
with the initial condition R(0) = 1, we have 

R(t) = (l-t/T ) 1 ^ a \ 
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Definition 2.1. Let fi G ^(M d ) be a radially symmetric probability 
measure. We define fx G ^([0, +oo)) by 

/x(J) = G M d : |x| G J}) 

for all Borel sets / in [0, +oo). 

In this section, we give a characterization of all the radially sym- 
metric first-kind similarity measure solutions to (11. ip in the case that 
c/ + a>5or<i-|-a = 4. We state the result as the following theorem, 
which generalizes Theorem 3.1 of [9] proved for the case a = 1. 

Theorem 2.2 (Characterization of similarity solutions). Let d > 3 

and VK(x) = x\x\ a ~ 2 . Assume that either a G [5 — d, 2) (and thus 
d > 4) or a = 4 — d. Then any radially symmetric first-kind similarity 
measure solution is of the form 



x 



R(t) d ^ \R(t) 
where 

fi = m 5 + m 1 5 pi 
for some constants mo, mi > and p\ > 0. 

Before we give a proof of the theorem, first we recall that for any 
radially symmetric measure solution fi t and t > 0, 

r°° x 
v(t,x) = - / ( j ) (\x\/p)p a - 1 dil t (p)— , (2.2) 
Jo M 

where <p : [0, oo) — > 1R is a function defined by 

0W= / , r ~ V ' da(y); (2.3) 

J St \ re i - y\ 

See, for example, [1] or [9]. From the dominated convergence theorem, 
it is easily seen that, when a G (2 — d, 2), the function is in fact 
continuous on [0, oo) and satisfies 

lim0(r)/r = O. (2.4) 

r— >oo 

Our proof of Theorem 12.21 is founded on the following key lemma. 

Lemma 2.3. i) Let d > 4, a G [5-d,2), andVK(x) = x\x\ a ~ 2 . Then 
the function <p(r)/r is a strictly decreasing function on (0, oo). 

ii) Let d > 3 ; a — 4 — d, and VK(x) = x\x\ a ~ 2 . Then the function 
<p(r)/r is a decreasing function on (0, oo). Moreover, it is a constant 
on (0, 1] and is strictly decreasing on [1, oo). 
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Remark 2.4. The proof of Theorem 14.11 of [9] uses the concavity of 
when a — 1; see Lemma 3.2 [9]. This property, however, is unavailable 
for general power a 6 (2 — d, 2) even under the condition d + a > 5 
or d + a = 4. Indeed, when a < 1 we have from (12. 3p that — >■ as 
r — )• oo. Since is strictly positive on (0, oo), it cannot be concave on 
(0, oo) in this case. 

Proof of Lemma \2.3[ We shall prove that under the conditions of the 
lemma, <fr/r is a C 1 function on (0, oo) and its derivative <fr'/r — <p/r 2 is 
negative on (0, oo). To this end, we rewrite ( 12.31) as 

r (r-cosfl)( sin0) rf - 2 



W oo d Jo A*-«{r,6) 
where 

A(r,6) = (1 + r 2 - 2rcos#) 1/2 

and Cc^ is the surface area of the unit sphere Si in M. d . For r e [0, 1) U 
(1, oo), a direct computation gives 



and 



r0'(r) — 0(r) 



(skitf)^ — (2 - a)r± '- A6 

" 1 (/i + / 2 ). (2.7) 



Integration by parts yields 



K ' A 2 ~ a J (d-l)A 2 ~ a 
d - 1 A 4 " Q v 



Clearly, 



I a = rj[ -(2-c)(sin^- 2 lr ™ a "' d6 
= r£-(2- a )(sin9r 2 A2 -J A f )2 a9. 



\2 
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r0'(r) — </>(r) 

_ r(2 - a)^ 
(d - l)u d 



d 



sin i 



(d-1) 



.sin i 



\<i-2 



A 2 " 



d0. 



(2.9) 



We shall first consider Assertion ii), then Assertion i). 
Assertion ii): d + a = 4. In this case, from (I2.9p we have 



r<f>'{r) - 0(r) = % °^ {d^ d - (d - l)* d _ 2 ), 

where 

9* = *k(r):=£^^M, k = 0,1, 2, ... (2.10) 
It is known that can be computed explicitly: 



W fc + 2 r„ nl 

on [0,1], 

r on (l,oo). 



Indeed, this is obviously true when k = 0. For k > 1, ujk+\^k{f) is 
the integral of \x — rei| _fc on the unit sphere in R fc+2 . Notice that for 
fixed r > 1, \x — re\ \~ k is a harmonic function on the unit ball in R fc+2 . 
Therefore, by the mean value theorem, we get Uk+\^k{ r ) — ^k+2 r ~ k - 
This also holds true when r = 1 by the dominated convergence theorem. 
The remaining case r e (0, 1) follows from the simple identity A(r, 9) = 
rA{l/r,6). Since Uk+2 = 2nuk/k, we get rcf>' — = on [0,1) and 
rcf)' —cf) < on (1, oo). Note that because \r— cos#| < A and | sin#| < A, 
the integral on the right-hand side of (12. 7p is absolutely convergent and 
continuous at r = 1 by the dominated convergence theorem. Thus (12.71) 
holds on the whole region [0, oo). So we conclude rcf)' — <f> = on [0, 1] 
and rcf)' — <f> < on (1, oo). 



H. DONG 



Assertion i): d + a > 5. In this case, we estimate I2 in a different 
way. For r G [0, 1) U (1, 00), integration by parts gives 

n Ft ■ „srf 2 r 2 -2rcosfl + (cosfl) 2 ,„ 
j 2 = _ r(2 _ «) J ( sm 6) d ~ 2 -p-^ 

= -r(2-a) j( (sin + ^ _ ^ J a (-2r + cos 6>) dfl 



-r(2-a) / (sin^ 2 -£- + 



A A ~ a (d-l)A A - a 
{4-a)r{sm6) d . o „. ,„ 

+ (d-l)W + 
Integrating by parts again in the last term above, we get 

% F, • md 2 ^ (*™6) d 
I 2 = -r(2-a) I (sm0) — z~ + 



A 4 " Q (d - l)A 4 " a 
2(4-a)r 2 (sing) d (4 - a)r ((sin 6) d+1 )' 
(d-l)A e ~ a + (d- l)(d+ l)A 6 -« 

x r, • mrf 2 r2 (sin#) d 
-r(2-a) / (sm^- 2 - — + 



A 4 ~ a (d - l)A 4 ~ a 
2(4 - a)r 2 (sin#) d (4 - a) (6 - a)r 2 (sin #) d+2 



d0. 



(c/-l)A 6 - a (d- l)(d + l)A 8 ~ a 

Therefore, by using (J27TJ) , (12. 8p . and the equality above, we obtain 



„ A Ud-i F ( sin e ) d ~ 2 (1 2 ( 4 ~ a ) ( sin 1 



r<//(r) - 0(r) = -H(2 - a)^ 1 



+ (4 r J a)(6 w" J ° ) yV- (2-11) 



(d- l)(d + l)A 4 

Since d + a > 5 and a < 2, it holds that 

(4-q)(d+l) 
~ (d- l)(6-a)' 

which yields 

2(4-a)(sin#) 2 (4 - a)(6 - a) (sin fl) 4 



1 



(d-l)A 2 (d-l)(d+l)A 

A- a /d + 1 2(sin#) 2 (6 - a) (sin #) 4 \ 
- d- 1 ' \6-a ~A 2 + (d+l)A 4 ) ~ °' 

The last inequality above is in fact strict almost everywhere because 
A > I sin 9 1 except at a single point where cos 8 = r. Therefore, we get 
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from ( 12~TTT) that r<j)'{r) - 0(r) < on (0, 1) U (1, oo). As before, the 
integral on the right-hand side of (12.71) is absolutely convergent and 
continuous at r = 1 by the dominated convergence theorem. Further- 
more, it is easy to check that 

lim{r0'(r) — <f)(r)} < 0; 

See also (I3.2p below. We thus conclude that r<p'{r)— 0(r) < on (0, oo). 
The lemma is proved □ 

We are now in the position to complete the proof of Theorem 12.21 

Proof of Theorem \2.2l With Lemma 12.31 in hand, the proof essentially 
follows an idea in the proof of Theorem 3.1 [9], by observing that, 
instead of the concavity of 0, the proof of Theorem 3.1 [9] only uses 
the monotonicity property of 0/r. 

We use a contradiction argument. Suppose (/if) is a radially sym- 
metric first-kind similarity measure solution such that there are two 
numbers 

< r\ < r 2 , rr, r 2 G supp p . 
Denote w(r) := |wo( re i)|- By 



w(r k )= / <j>(r k /p)p a - 1 dfk(j>), A; = 1,2. 
J o 

It follows from Lemma 12.31 that 

-0(n/p) > -<t>(r 2 /p) Vp G (0, oo), 
n. r 2 

with the inequality being strict for p in a small neighborhood of r\ 
because r 2 /ri > 1. Since r\ G supp/io, we get 

— iw(ri) > — w(r 2 ). 
n r 2 

However, from ( 12. ID w is a linear function on supp (iq: 

— w(ri) = — w(r 2 ). 
r\ r 2 

Therefore, we reach a contradiction. The theorem is proved. □ 
3. The case d + a e (4,5) 



In Theorem 12.21 we proved that double or multi delta- ring similarity 
solutions to the aggregation equation (11.11) solution does not exist when 
a+d> 5 or a+d = 4 by showing that any radially symmetric first-kind 
similarity solution is of the form 

Po = m 5 + mi^. 
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In this section we show that this characterization result is still true in 
the case a + d G (4, 5). 

Theorem 3.1. Let d > 3 be an integer, a G [4 — d, 2), and VK(x) = 
x\x\ a ~ 2 . Then any radially symmetric first-kind similarity measure so- 
lution is of the form 



for some constants m ,mi > and p\ > 0. 

In fact, we shall prove that the claim of Lemma 12.31 i) still holds 
if we only assume d > 3 and d + a > 4. As soon as this is verified, 
one immediately proves Theorem 13.11 by using the very same proof of 
Theorem 12.21 It seems to us that the remaining case d + a G (4, 5) is 
more involved. One of the difficulties to adapt the proof of Lemma 12.31 
is as follows. When d + a > 5, the strict negativity of r<p' — is proved 
by obtaining a complete square in the integrand. This method does 
not work when d + a G (4, 5), since in the limiting case d + a = 4 the 
function r<j)' — is equivalent to zero for r G [0, 1]. 

Here we prove Lemma [2.31 i) for d > 3 and d + a G (4, 5) by directly 
evaluating the coefficients of the series expansion of (12. 9p . In the sequel, 
we always assume d > 3 and d + a G (4, 5), and we write d + a = 4 — 7 
for some 7 G (—1, 0). 

First of all, we observe that, for any r G (0, 00) and 9 G (0, 7r), 




where 



m 5 + mi8 



in 



A(r,9) = rA(l/r,0). 



Thus, we have 




and 




Therefore, by ( 12. 9j) , to show 




•(r0'-0)<O in [0,oo), 



it suffices to prove 



^(r)<0 Vre(0,l]. 



(3.1) 
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We shall prove ( 13. ip in the remaining part of the section. Clearly, by 
using the property of fy k defined in (I2.10p . we have ?p(0) = 0. At r = 1, 
we have A(r,9) = 2sin(#/2) for 9 G [0,7r]. Thus, we compute 

r/ (sing)" (sing)^ 2 \ AQ 



2 ^r( dcosd ^- (rf - i)cos "©) si 

2^(dB( d 4±, 1 -^-)-(d-l)B' d - 1 l ~ 1 



[ d -\ 



2 ' 2 

7(rf - 1) Bf^,l^)<0. (3.2) 



<i — 7 

For r G (0, 1), we use change of variables to get 



fsin#) d .fsin^- 2 



rf (l_ t 2 )( d-l)/2 l)(l_ t 2 )( d-3)/2 



J_ x \{\ + r 2 _ 2rt)( rf +7)/2 (1 + r 2 - 2rt)( d +^ 2 )/ 2 

1 /-(l+r) 2 . 

2r 



df 



rfs -( d+7 )/2^_ (1+r 2_ s) 2 /(2r) 2y 

(l-r)2 ^ 

- (d- l) S -^" 2 )/ 2 (l - (1 + r 2 - S ) 2 /(2r) 2 ) { "" 3)/2 ) ds. 

The last integral can be evaluated explicitly. We expand the outcome in 
a series with respect to r. After a direct but quite lengthy computation, 
we obtain that, for d > 3 odd, 

M = -*± [i+ w$-w ^ (3 ' 3) 



and for d > 4 even, 
where 

ai = - 7 (2- 7 ), a k+1 = a k (2k + ^(d + 2k -2 + 7), fc = l,2 



^(r) = -7T V rr:a k r 2 \ (3.4) 

^ (rf+2A;)!!(2A;-2)!! ' v ; 



It is clear that, for 7 G (—1,0), the series above are absolutely conver- 
gent in [0, 1) and all the coefficients are strictly negative. Therefore, 
we conclude ip{r) < in (0, 1) and Theorem 13.11 is proved. 
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Remark 3.2. We note that the argument above does not apply to the 
case 7 < — 1 . Thus it seems unlikely to unify the proofs of Lemma 12.31 
for different regimes of a in this section and Section [2j 

4. Existence of multi delta-ring solutions 

In this section, we shall construct multi delta-ring similarity solutions 
to (11. ip when d + a < 4. Let n > 2 be an integer. The n delta-ring 
solutions are of the form 

/ x 



R(t) d ^ \R(t) 
where 

n 
k=l 

for some positive constants rrik, pk, k — 1, . . . , n satisfying 

n 

y^m fc = 1, < pi < p 2 < ■ ■ ■ < p n < oo. 

k=l 

In the 2D case with a = 1, the existence of two delta-ring solutions 
was obtained in [1] by using a fixed point argument. The main result 
of this section is the following theorem, which gives the existence of n 
delta-ring solutions for any n > 2 when a + d < 4. 

Theorem 4.1 (Existence of multi delta-ring similarity solutions). Let 

d > 1 be an integer and a < 2 such that a G (2 — d, 4 — d) . Then for 
any integer n>2, there exists an n delta-ring solutions to (II. ip . 

The lemma below will be used in the proof of Theorem 14.11 We 
postpone the proof of it to the end of this section. 

Lemma 4.2. Let d > 1 be an integer and a < 2 such that a G (2 — 
d, 4 — d). Then we have 

0(1) > 0'(O) > 0. (4.1) 



Proof of Theorem 4-1 We shall construct a certain type of delta-ring 
solution where the radii of rings of the initial data are powers of a given 
large number. Our goal is to determine how to distribute the unit mass 
on these rings, so that the solution will be a similarity solution of the 
first kind. From the proof below, it should be apparent that for the 
existence of such a distribution, all we need is that ratios of different 
radii are sufficiently large (or small). Of course, there could be many 
other multi delta-ring solutions which are not of this structure. 
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Let A be a sufficiently large constant to be specified later. We take 

Pk = A fc_1 , m k = p 2 k ~ a rh k , k = 1, . . . , n. 

Again, we denote w(p) = \v(pei)\. Thanks to (12. ip . in this case a 
necessary and sufficient condition for p t to be a similarity solution is 
that 

w(p k )/p k = w(p j )/p j , kj = l, ...,n. (4.2) 

So to prove the theorem it suffices to show that, for sufficiently large A, 
there exists fh k , k — 1, . . . , n such that the condition (14. 2ft is satisfied. 
From (12.21) . we have 



w{Pk) = ^2<P(Pk/Pj)p* X m i = J^0(A fc ^Pjrhj. 
3=1 j=i 

Therefore, 

n 

w(p k )/p k = ^(j){\ k - j )\ j - k rhj, 

3=1 

and ( 14.21) is equivalent to the existence of a solution (xi, . . . ,x n ) to the 
linear algebraic equation 

fi A (x l7 ...,x n ) T = (l,...,l) T (4.3) 

satisfying Xj > 0,j = l,...,n, where B\ = [Bx,kj]jk=i is an n x n 
matrix and 

Bx, kj = 0(A fc - J )A J - fc 
It follows from O that 

lim Bx,kj — -Soo.fej := 

A— >oo 

which together with Lemma 14.21 implies that is upper-triangular 
and non-degenerate. It is a simple fact that the equation 

B 00 {x l ,...,x n ) T = (1,...,1) T 

has a unique solution which is given by 

( Xl , . . . ,x n f = (0(1))~ 1 ((1 - cT~\ (1 - c)- 2 , . . . , 1) T , 



3,k = 


1,.. 


.,n. 


m 


if 


k<3, 


0(1) 


if 


k = j, 





if 


k > j, 



where c = <j>' (0) / <j>(l) . By Lemma 14.21 we have 1 — c > 0. Since 
.Boo is non-degenerate, by a standard continuity argument, for A > 
sufficiently large, (14.31) has a unique solution, which satisfies Xj > 0, j = 
1, . . . , n. The theorem is thus proved. □ 

We finish this section by giving the proof of Lemma 14.21 
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Proof of Lemma \4-2\ First we remark that (14. ip essentially follows from 
the calculations in Section [3J Indeed, by the expansions (13.3j) and (13.4)) . 
one can easily see that ip > in (0,1), which implies that 0(r)/r is 
strictly increasing in (0, 1]. Below we give an alternative proof without 
using the complicated series expansions. 

As before, denote 7 = 4 — d — a G (0, 2). By (12. 6p . we have 

0'(O) = r(sin9) d - 2 (l - (2 - a)(cos#) 2 ) d6 

Jo 



UJd 



B( d -^±A- ( 2-o)B< d - 1 3 



2 '2/ v y V 2 ' 2 



w d _i2-7^/d-l 1 



u d d V 2 2 



w d d V 2 / V2/V V2 
where B(-,-) is the Beta function. Here, we used the formula 

^(cosensmeyae = ^(£±1 v /?, 7 > -1. 

By (12. 5p . we compute 

Ud-i f T {l-cos6){sm0) d - 2 ^ 



u d J (2 - 2cos#) 1 - Q / 2 
Jo 2 



; (2 sin(0/2) cos(fl/2)) d 2 (2 sin(fl/2)) a dfl 



-2 1 " 7 (sin(0/2)) 2 ^(cos(0/2)) < *~ 2 d0 
Jo 



Ud 

w d -i_ 2 i_ 7jB ^3-7 d-1 



w d V 2 2 

^-rfi^W^ fr^-^ 2 ))' 1 . (4.5) 



It follows from the duplication formula for the Gamma function that 
'3-7\ m ,„/lW. \M2-T 



2-T(i)r(2- 7 )(r(" 



2 

This together with (14.51) gives 



(4.6) 
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By using (14.4p and (14.6p . to prove ( 14. ip it suffices to show 
which is equivalent to 
Since 

+ 2 rf-7 + 2 4-7 
2-7 + ^— = i + 



and 



(2 _ 7) .l±l<< i -^ 2 



2 2 2 ' 

the inequality (14. 7p is an immediate consequence of the log-convexity 
property of the Gamma function. □ 

5. An EXAMPLE 

Finally, we present an example of first-kind similarity solutions with 
support on open sets in the limiting case a = 2 — d, which generalizes 
a ID example constructed in [1] with a — 1. 

Example 5.1. Assume d > 2 and a — 2 — d. Let Mo(x) = rfcjj 1 /^^!. 
Then u given by ( II. 2p with R{t) = (1 — tud) 1 ^ is a measure similarity 
solution to the aggregation equation ( 11. ip . Moreover, T = c^ 1 is the 
blowup time of the solution. 

As mentioned before, to prove that u is a similarity solution, it suf- 
fices to show that 

poo 

w (r) ■= / <j)(r/p)p a ~ 1 dfio(p) 
Jo 

is a linear function with respect to r on [0, 1), which is the support of 
/t - For a = 2 — d, from (12. 5p we have 

0(r) ~/o AW) ^ 

which can be evaluated explicitly: 

on [0,1), 

0(r) ={1/2 at r = 1, 

- 1 - d on (l,oo). 
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Consequently, for r G [0, 1), we have 

w(r) := [ 1 <f>(r/p)p a - 1 dfL (p) 
Jo 

= f <j ) (r/p)p a - 1 dfi (p) 
Jo 

= r^ d p (Br) = ^r. 

Therefore, indeed w is a linear function on [0, 1) and a is a similarity 
solution to (11. ip . Furthermore, we have A = Ud/d. Following the 
argument at the beginning of Section [2} we obtain 

T = (2 - a)" 1 A" 1 = ooj\ R(t) = (1 - tu d ) 1,d . 
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